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Dark matter (DM) without a matter asymmetry is studied in the context of Twin Higgs (TH) 
theories in which the LHC naturalness problem is addressed. These possess a twin sector related to 
the Standard Model (SM) by a (broken) Z 2 symmetry, and interacting with the SM via a specific 
Higgs portal. We focus on the minimal realisation of the TH mechanism, the Fraternal Twin Higgs, 
with only a single generation of twin quarks and leptons, and SU (3)' x SU (2)' gauge group. We show 
that a variety of natural twin-WIMP DM candidates are present (directly linked to the weak scale 
by naturalness), the simplest and most attractive being the r' lepton with a mass m^i > mniggs/S, 
although spin-1 W'^ DM and multicomponent DM are also possible (twin baryons are strongly 
disfavoured by tuning). We consider in detail the dynamics of the possibly (meta)stable glueballs in 
the twin sector, the nature of the twin QCD phase transition, and possible new contributions to the 
number of relativistic degrees of freedom AA^eff. Direct detection signals are below current bounds 
but accessible in near future experiments. Indirect detection phenomenology is rich and requires 
detailed studies of twin hadronization and fragmentation to twin glueballs and quarkonia and their 
subsequent decay to SM, and possible light twin sector states. 


I. INTRODUCTION 

Models based on the Twin Higgs (TH) mechanism 
address the LHG fine-tuning problem and solve 
the little hierarchy problem by introducing a twin sector, 
that is, in its simplest realisation, a copy of the Stan¬ 
dard Model (SM), regarding both field content and in¬ 
teractions. At tree level, the Higgs sector of the theory 
respects a global SU{A) symmetry (in fact an 0(8) sym¬ 
metry in the most attractive cases m) acting on the com¬ 
ponents of the pair of Higgs doublets where H' 

is the twin Higgs doublet (throughout primes denote ob¬ 
jects in the twin sector). A discrete Z2 between the two 
sectors ensures equality of their couplings, which results 
in ^^(dj-symmetric radiative corrections to the Higgs 
mass squared. The SU{A) symmetry is broken at one 
loop order by radiative corrections to the Higgs quar- 
tic coupling and the SM Higgs is realised as a naturally 
light pseudo-Nambu-Goldstone boson of the approximate 
SU{A). The Z2 symmetry needs to be broken, explic¬ 
itly or otherwise, for the SM Higgs to acquire a phe¬ 
nomenologically viable vacuum expectation value (vev), 
for an exact Z2 would imply that the vev’s in the two 
sectors are equal, a possibility that is excluded by Higgs- 
coupling measurements. Denoting the SM and TH vev’s 
as u « 246 GeV and / respectively, the fine-tuning aris¬ 
ing from this difference of vev’s is ^ 2v^/p, i.e. a mild 
^ 20% tuning for the minimum experimentally allowed 
ratio //u « 3. The physical light Higgs state, h, is shared 
between the SM and twin sectors with couplings to SM 
and twin sector states modified by cos{v/f) ~ l — v'^j2p 
and siii{v/f) ~ v/f respectively. It is important to also 
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bear in mind that the TH theory needs to be UV com¬ 
pleted at some cutoff scale Muv < 47r/ (for definiteness 
we take Muv ~ 5TeV). 

For the TH mechanism to operate, the twin sector does 
not need to be an exact copy of the SM, a reduced field 
content sufficing. This simplified version - the Frater¬ 
nal Twin Higgs (FTH) [5] - has as its minimal ingredi¬ 
ents twin weak SU{2y and colour SU{3y interactions, 
and a twin third generation consisting of top and bottom 
quarks Q', and a lepton doublet L' required by 
SUiyiy anomaly cancellation, and a twin Higgs doublet 
H'. Right handed twin leptons may be added to the the¬ 
ory rendering the leptons massive, although they are not 
required by the TH mechanism. For the TH mechanism 
to still be effective without introducing significant extra 
tuning, the twin top Yukawa yn can only differ by at 
most 1% from yt for Muv ~ 5 TeV. The need for gauged 
twin SU{3y then becomes apparent: radiative correc¬ 
tions to yt and yp could make them differ significantly at 
the weak scale even if they coincided at Muv- A gauge 
coupling ^3 differing by less than ~ 15% from 53 at the 
cutoff ensures that the running of yt and yp is close 
enough. However, the running of the 53 and (73 couplings, 
and thus the dynamical scales, differ because of the dif¬ 
ferent field content and masses of the two sectors. With 
only one quark generation in the twin sector, and for the 
allowed range of values of g '^, we have a twin confinement 
scale A(;jcd ~ 0-5 “ 20 GeV [6]. (Unless otherwise stated 
we will take our default value as Aqj^^q « 3 GeV; the 
preference for this choice, or larger, is justified in Sec¬ 
tion VHI.) The presence of an 5't/(2)' gauge group in the 
twin sector also follows from fine-tuning considerations, 
with the twin <72 coupling allowed to differ by at most 
~ 10% from the SM value. A gauged 1/(1)' is not neces¬ 
sary from a naturalness perspective, although it remains 
an accidental global symmetry of the twin sector. 

In this paper, we explore the possibilities for dark mat¬ 
ter (DM) in TH scenarios, focusing on the minimal FTH 
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models in the case where a matter-antimatter asymme¬ 
try in the twin sector is not present and where U{iy is 
a global symmetry. (We reserve the study of both twin 
asymmetric DM, and the effects of a gauged U{iy, for 
a companion paper [7].) Crucially for our later discus¬ 
sion, intrinsic to the success of the TH mechanism is the 
fact that the twin sector interacts with the SM via the 
Higgs portal with a strength determined by f /v, leading to 
specific predictions for DM signals. In addition, we will 
show that the most attractive twin DM candidates in the 
absence of an asymmetry are the twin leptons, which nat¬ 
urally have a ‘twin-WIMP miracle’ as they freeze out via 
twin weak interactions whose strength is set by — 92 
and G'p = {v/f^Gp both of which are directly tied to 
SM weak interaction values by naturalness. Since the 
Higgs portal forces the SM and twin sector to be in equi¬ 
librium until temperatures well below the EFT cutoff, 
we have a purely thermal freeze-out scenario, at least for 
states that are not UV relics. 

As asymmetric DM is not our concern here, we will 
primarily focus on the heavy twin quark limit, i.e. 
mt ', mp 3> Aq(^q , which arises naturally for yp = yt 
and yp « ys and for the values of f /v > 3 that are al¬ 
lowed. Collider signals of the TH model in this regime 
were treated in detail in [5], and also considered in [5]. 

Finally we remark that the idea of a mirror world, ei¬ 
ther exact or partial, has a long and rich history [9Hn], 
see e.g. [T^]. (For a recent review of aspects of mirror 
world physics we refer the reader to [13] .) Often such the¬ 
ories lead to a variety of interesting DM candidates with 
overlap with those studied here. Here we are considering 
a particular, approximate mirror world interacting with 
the SM via the Higgs portal, as directly motivated by 
the LHC naturalness problem. Previous investigations 
of symmetric dark matter candidates in models with a 
similar philosophy include mi and ca¬ 


ll. STABLE & METASTABLE TWINS 

At temperatures well below the SU{2y phase transi¬ 
tion, where anomaly effects are exponentially small, the 
FTH model has an accidental C/(l) global symmetry as¬ 
sociated to twin baryon number Bb If Majorana mass 
terms for the right-handed twin leptons are forbidden 
(for example, due to a discrete symmetry), then there 
are also accidental C/(l) twin lepton number and twin 
‘charge’ symmetries, with associated conserved numbers 
L' and Q'. Ultimately, we might expect these global sym¬ 
metries to be explicitly broken by higher dimensional op¬ 
erators, possibly connected with Planck-scale physics, or 
by terms from the UV completion of the TH models that 
connect the SM and twin sectors in ways beyond the TH- 
mandated Higgs portal interaction itself. For the pur¬ 
poses of this work we assume that these new interactions 
are sufficiently weak that the lightest states carrying B', 
L' and Q' are stable on timescales ^ 10 ^although 
decaying DM is a natural possibility in TH models. 


The discrete symmetries P and C in the twin sector are 
maximally violated by 311(2)' interactions but, in princi¬ 
ple, GP can remain conserved. Although we focus on the 
CP-preserving scenario, breaking of GP in the twin sec¬ 
tor, due, say, to an un-cancelled 0 qcd' term, is allowed 
and can have important consequences. We mention the 
changes this makes where appropriate. 

We first consider the case where < to,-', 

m ^>, mr> < mw ', implying that both v' and r' are 
stable states and therefore automatic DM candidates. 
We will focus on the regime of heavy r' (m,-' ^ m/i/2), 
whereas we will allow v' to be heavy or effectively mass¬ 
less. When v' is effectively massless it will behave as 
dark radiation (DR), contributing to the number of ef¬ 
fective neutrinos, Aeff. (We discuss the issue of Aes in 
detail in Section VHI.) If -I- m,-' > rn\y'-, the twin 
W'^ gauge bosons will also be stable, as automatically 
nit' > mw' in TH models, and therefore could 
contribute significantly to the DM density. 

In the strong sector, things are more involved. One 
stable state is the lightest twin baryon. A', made out of 
three b' quarks in a spin 3/2 state. In the absence of an 
asymmetry in the twin sector, b'b' pairs annihilate effi¬ 
ciently into gluons rendering A' irrelevant as a DM can¬ 
didate unless mb' ^ 1 TeV, a case we discuss in Section V 
(if the freeze-out temperature of the b' quark-antiquark 
annihilations is below the phase transition temperature 
then A' (anti)baryons annihilate efficiently into glueballs 
and quarkonia). However, twin QCD glueball states are 
themselves of potential interest depending upon the UV 
completion of the TH theory. In the heavy quark regime 
the spectra of glueball and quarkonia are well separated 
(and relatively well known). The lightest glueball, whose 
mass is mo « G.SAqqq (THldT], has quantum num¬ 
bers 0'^'*" and therefore mixes with the Higgs, with mixing 
angle |3] 

^ a'^vTo _ vml 
67r/2(m^ —mg) 

with To the 0++ glueball decay constant given by Pq = 
3.06mo/(47ra3) [T7| and we have assumed m^ » m^ in 
the final step. It therefore promptly decays to light SM 
states (to++ ^ 4 x 10“^° s if we take Aq(2J3 = 3 GeV and 
//w = 3). In the case where the twin neutrino is light 
compared to Aq(^q, all other glueballs decay, via 311(2)' 
and heavy quark-loop induced interactions, to (eventu¬ 
ally) combinations of the O’*""'' glueballs and GP in ap¬ 
propriate angular momentum states, so leaving no stable 
twin-QCD states apart from A'. 

In the case where both t' and P are heavy, such 
that twin-lepton pairs are not kinematically accessible 
in glueball decays, another two glueballs become po¬ 
tentially relevant: a 0 ^ glueball, with mass mQ-+ « 
l.Smg [T6l[T7], and a 1"' glueball, with mass mi+- « 
1.7mo USKIT]. These glueballs can be potentially long- 
lived metastable states, possibly with cosmologically long 
lifetimes, though this sensitively depends on the issues 
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of additional SM-twin-sector interactions arising from 
the UV completion and/or twin CP-violation. As we 
discuss in Section VI the freeze-out abundance of these 
glueballs is very small, ~ 10“®pdm, so they are not sig¬ 
nificant gravitationally, though they may be constrained 
by CMB or cosmic ray observations if their lifetimes are 
long. Also, as we discuss in Section VII, both 0 ^ and 
I^ glueballs have the potential to produce novel indirect 
detection signals if they decay with lifetimes in the ranges 
T ~ 1 —10^ s or T ^ 10^^ —10^^ s. Whether such lifetimes 
are achieved depends on the value of Aq^^q, and on the 
UV completion, a subject we reserve for Section VI. 


III. TWIN SU{3) PHASE TRANSITION 


Before we can calculate the relic densities of the stable 
twin states, we first must consider the nature of the twin 
QCD phase transition, and whether it leads to significant 
dilution of relics by entropy production. 

If Wb' ^ 8Aq(^q , then the one-dynamical-quark-flavour 
QCD' phase transition is a smooth cross-over [T8ll20j . 
with no significant non-equilibrium dynamics. As ruf,/ is 
increased, the transition becomes second order at a crit¬ 
ical value ^ SAqcd' (it should be noted that this upper 
limit is potentially uncertain by a few times Aqj^q), and 
above that is (weakly) first order, as demonstrated by 
analytical arguments [niEi] and lattice studies [22]. 

Investigating the dynamics of the first-order phase 
transition in the heavy-quark case, lattice studies give 
a critical temperature of ~ I.26Aqqq, |53|, bub¬ 
ble wall surface tension a ~ O.OISST/*, and latent heat 
Pl — 1-4T/ [231 . Since the pressure is continuous across 
the phase transition, the pressures of the confined and un¬ 
confined phases are equal at Tc- If the unconfined phase 
manages to supercool to a temperature T = (1 — d)Pc, 
then a bubble of confined phase can grow due to the pres¬ 
sure difference between the phases, if it is large enough 
to overcome the surface tension, i.e. if it has radius 
R> Rc = The free energy cost of a ‘critical bubble’ 
of radius Rc is 

AP, = AnRla - ^ttRIAP = ^ . (2) 

Assuming that the supercooling is small, i5 <C I, so that 


AP = Pg{T) - Pg{T) ~ -STaisoin) - SgiT,)) = 5 pl . 

( 3 ) 

where Pg is the pressure in the confined phase, Pg is 
the pressure in the unconfined phase, and SG,Sg are the 
corresponding entropy densities, we have 


AFn IGtt 


T 


3 piTc 


~ 3 X 


-5r-2 


( 4 ) 


The rate per unit volume of bubble nucleation due 

to thermal fluctuations is set by T/V ^ 

pS] . This becomes comparable to i.e. one event 


per Hubble volume per Hubble time, when AFc/T < 
log{T^/F[^). Assuming that d ^ 1, this is equivalent 
to 3 X 10-5d-2 < AlogiTc/H{Tc)) ~ 160, so S > Sn = 
4 X 10“"^. Since grows by an e-fold with a ~ 10“® 

drop in log T, the nucleation rate quickly becomes large 
as T drops. 

This extremely large nucleation rate means that, once 
6 is somewhat larger than d„, the latent heat released 
from bubble nucleation and expansion must heat us back 
above the nucleation temperature, or the transition must 
complete entirely in a very small fraction of a Hubble 
time (in the most extreme limit, for 6 ~ 10(5„ we have 
g-AF^/T ^ 0(1), and nucleation effectively happens ev¬ 
erywhere at once). The rapid-completion case would 
arise if e.g. the heat capacity of the supercooled uncon¬ 
fined phase was very high, so that the expansion neces¬ 
sary to cool it from Tc to (1 — 5)Tc removed almost all of 
the latent heat of the transition. In the former case, the 
confined phase bubbles grow, releasing latent heat, until 
the temperature has been brought back up to where 
the pressures are equal. Once there, the confined phase 
bubbles can only grow as Hubble expansion removes en¬ 
ergy from the mixture. In this way they grow until they 
occupy the entire volume, with the temperature remain¬ 
ing constant at T^. 

In either case, the out-of-equilibrium part of the phase 
transition, during which there is a pressure difference 
driving confined phase bubble expansion, occurs at tem¬ 
peratures only very slightly below T^, so can only pro¬ 
duce a small amount of entropy. The maximum possi¬ 
ble increase in the overall entropy density occurs in the 
rapid-transition case, where the pressure-driven growth 
transforms all of the volume from the unconfined to the 
confined phase, in which case the overall increase in en¬ 
tropy density is 

As ~ ^ ~ few X ~ I0"^r/* . (5) 

In the (naively more likely) case where most of the tran¬ 
sition occurs in quasi-equilibrium, the maximum entropy 
production is even smaller, As cx S^p^lTc ^ 2 x I0“^T^, 
since only a fraction oc of the volume is converted by 
out of equilibrium growth. Either way, there is no signif¬ 
icant effect on relic densities. If the transition is simply 
a cross-over, as occurs for light enough 6', then there is 
again no significant entropy production. For a crossover 
or a quasi-equilibrium transition, there will not be sig¬ 
nificant gravitational radiation production, rendering the 
FTH model less encouraging for gravity wave signals than 
the other scenarios discussed in |26j . 

We also remark that if twin-CP is violated by dqcD' 

0 in the pure SU (3) case with no light quarks, the tran¬ 
sition remains first order but with parameters such as 
the critical temperature, Tc{0), and the latent heat, now 
depending on 0 qcd' in a periodic fashion but otherwise 
poorly constrained, either analytically m or from lattice 
studies [55] ■ These studies indicate that for dqcD' ‘C tt 
the transition temperature decreases Tc{0) < T'c(O), while 
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the strength of the transition slowly increases. However 
the properties of the phase transition at 0 qcd' 7 ^ 0 are 
not firmly established, especially once 0 qcd' is not small. 
Thus to be conservative, in the cases with twin-CP vi¬ 
olation we take 0 qcd' tt, a limit that is sufficient for 
our purposes. 

IV. TWIN TAU DARK MATTER 

Having argued in Section HI that the twin-QCD phase 
transition leads to no significant dilution of relics by en¬ 
tropy production we now proceed to calculate the freeze- 
out density of the stable twin-sector states. We start with 
the simplest case, that of the twin t lepton (since 17(1)' 
is not gauged, the situation for v' is identical to that of 
t' assuming a suitable Yuakawa coupling giving a Dirac 
mass). In most of parameter space the annihilation of 
t'’s dominantly proceeds via twin-S'{7(2) weak interac¬ 
tions into the (assumed lighter) &'-quarks/quarkonia and 
v'v' pairs. Annihilation via the Higgs, with couplings 
that are given by is subdominant apart from 

a narrow resonance region around vrir' ~ mh/‘i- 

Figure]^ shows the contribution to the present energy 
density of the Universe from r' species, normalized to the 
observed DM density for different values of f /v. This 
calculation, along with the other relic density calcula¬ 
tions throughout this paper, was performed using the 
MicrOMEGAs software package [5S]. For f /v ~ 3, the 
least tuned case, the observed DM density is obtained 
for rur' ~ 63GeV. Larger values of the ratio f /v, which 
imply worse tuning, result in larger DM massesjj 

We emphasise that, for (symmetric) twin DM candi¬ 
dates with 0(0.1) Yukawa couplings, we can obtain the 
correct relic density since the couplings and mass of the 
W bosons are set by the TH mechanism to be §2 — 52 
and Mw' — {f/v)M’^±. These are constrained by natu¬ 
ralness to be close to SM weak interaction values, giving 
rise to a natural ‘twin-WIMP-miracle’. 

Turning now to direct detection, scattering of r' with 
SM nuclei occurs, at tree level, via Higgs exchange, and 
this process sets the scattering cross section in direct de¬ 
tection experiments. In Figure]^ we show the spin inde¬ 
pendent scattering cross section per nucleon off SM nu¬ 
clei for t' dm as a function of and for values of f /v 
such that the correct DM abundance is obtained. For 
f/v ~ 2> — b (tuning 20—8%), the predicted direct detec¬ 
tion signatures fall below current experimental bounds 


^ For concreteness, we take my Ri 15GeV, which saturates the 
experimental bound coming from constraints on the Higgs width 
for f/v ^ 3. This also corresponds to the Z 2 symmetric value 
(f/v)mf, for f/v ^ 3. The twin bottom Yukawa in the FTH 
is only constrained by tuning to be <g yy. Also, as long as 
my , Aqqq <IC m^i , the t' relic density is insensitive to the pre¬ 
cise b' mass. Thus, different (but still sufficiently light) values of 
my would not affect our conclusions. 



FIG. 1. Contribution to the energy density of the Universe 
from t' species normalized to the observed DM energy density 
as a function of m.^' for different values of f/v. Light (dark) 
pink area indicates the 2-sigma bounds from invisible Higgs 
width and modified couplings to visible sector particles in the 
case f/v = 3 (3.5), whereas f/v = 4,5 remain unconstrained 
in the region of parameter space shown. Note that, if 
is large enough so that mo > 2my , then annihilations of low- 
mass t's will have significant non-perturbative corrections. 
However, this regime generically leads gives too high a t' 
density, so is not of primary concern here. 

[50] but above the region of parameter space that will be 
probed in the very near future by LUX |31j . Larger (more 
tuned) values of f/v will be probed by next-generation 
experiments such as LUX-ZEPLIN (LZ) [32] . 

V. MULTICOMPONENT, W, & A' DARK 
MATTER 

In the case where the sum of r' and v' masses is larger 
than the W mass, the latter is not able to decay. In 
the regime where to,-' ^ m^' and m,-/, m^i < mw ', this 
implies that all three states are stable and may signifi¬ 
cantly contribute to the DM energy density, opening a 
possibility for a 3-component DM scenario. 

Figure[^shows the contribution to the DM energy den¬ 
sity of these three particle species (normalized to the ob¬ 
served value) for different values of the twin weak cou¬ 
pling, that we allow to vary by 10% from its central 
value 52 = 52 ~ 0.64. For concreteness, we have taken 
rUr' = rn^i « 0.55 mw', with mw' = 52//2- As one can 
see, the observed DM energy density is only achieved for 
relatively large values of the ratio //n, where the fine- 
tuning is in the range 5% to 1%. This occurs since r' 
and v' are forced to be heavier than considered in the 
//u ~ 3 case, so their annihilation cross sections set by 
m^//^ are larger — to compensate, / must be increased. 

As can be read off from Figure]^ the correct DM abun¬ 
dance is obtained for f /v ~ 9.7 for 52 = 52 j which implies 
a tuning of approximately 2%. In this case, t' and v' 
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FIG. 2. Dashed green line is the spin independent scattering 
cross section per nucleon for r DM for those values of rriT-i 
and //w such that the correct DM density is obtained. Red 
dots point out particular values of //u (indicated in num¬ 
bers). Blue area: LUX current bounds |30| : blue line: LUX 
projected sensitivity (300 live-days) [3T]; pink area: region of 
parameter space ruled out by bounds on the invisible Higgs 
width, and on modified couplings to visible sector particles [6] . 


species would contribute roughly 75% to the DM energy 
density, with W' species making for the remaining 25%. 

Regarding direct detection experiments, the predicted 
spin independent scattering cross section per nucleon for 
all three particle species is of order ~ 10“^®cm^ for values 
of the masses that result in the correct DM abundance. 
This lies around an order of magnitude below LUX cur¬ 
rent projected sensitivity for the range of masses consid¬ 
ered, which means that next-generation direct detection 
experiments such as LZ [32] will be able to cover the 
relevant region of parameter space. 

Small variations of the values of t' and v' masses 
around the case we have considered do not make a 
significant difference to our conclusions, except when 
rrir' ~ rn^' but rrir' + ra^/ < mw'- In this case, the 
W is no longer stable and then only r' and v' species 
would contribute to the DM density. In this 2-component 
DM scenario, sufficient annihilation requires m,-' and 
rriv' to be in the mass range above mhl‘2-, automatically 
evading invisible Higgs width constraints. The differ¬ 
ent contribution to the DM density from the two par¬ 
ticle species would depend solely on the ratio of their 
masses: for equal masses, both components would con¬ 
tribute 50%, whereas if they differ by approximately 
10 GeV the right DM abundance requires rriu' ~ 70 GeV 
(therefore m,-/ « 80 GeV) and v' and r' species would 
make for 65% and 35% of DM respectively. Regard¬ 



FIG. 3. Gontribution to the energy density of the Universe 
from r', u' and W', species normalized to the observed DM 
energy density, as a function of f/v and for different values 
of the twin weak coupling g 2 . Vertical lines represent tuning 
contours. We indicate the W mass for the three different 
values of considered when the right DM density is achieved. 


ing direct detection signals, this 2-component scenario 
is analogous to the single-component case discussed in 
Section IV, for interactions between the DM species and 
the visible sector proceed only via Higgs exchange. 

Finally we turn to the most complicated of the possible 
twin DM candidates, the A' baryon. Although for light h' 
quarks, and in the absence of a matter-antimatter asym¬ 
metry (a subject of a companion paper [ 7 ]), the spin-3/2 
A' baryons efficiently annihilate to glueballs and quarko- 
nia, leaving an uninterestingly small freeze-out density, 
this is no longer the case if the 6'-quarks, and thus the A' 
baryons, are sufficiently heavy, mb' ^ 1 TeV Aqqq. 

To estimate the freeze-out density of such states let 
us consider the case where the freeze-out temperature is 
well above Aq(^q, a situation that applies if mb> is suffi¬ 
ciently large. In this case we may self-consistently work 
with h' quarks and twin-gluons, and first calculate the 
freeze-out density of b' quarks, via a leading annihilation 
cross section to two gluons that parametrically goes as 
av ^ {a'^/mb')'^■ We find that a numerical evaluation of 
the annihilation rate leads to a freeze-out temperature 
Tf ^ mb '/30 and gives a substantial freeze-out density of 
6' quarks and anti-quarks only once mb' ^ 1 TeV, which 
implies very large f/v ^ 30 and thus a very badly tuned 
TH model (here we have taken yb' ^ 0.2 so as not to 
have yet further tuning at 1-loop). We therefore come 
to the conclusion that the A' baryon in the absence of a 
matter-antimatter asymmetry is a poor DM candidate in 
TH models. (We remark that the freeze-out density esti- 
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mated this way provides the most optimistic estimate of 
the minimum b' mass. The reason is that the 6' quark and 
anti-quark densities do not simply translate, via a factor 
of 1/3, into the final freeze-out density of A' baryons and 
anti-baryons. Post twin-confinement only a proportion 
of 6'-(anti-)quarks end up in A' (anti-)baryons, and thus 
are stable relics, compared to the number that form b'b' 
quarkonia and therefore quickly decay. We expect this 
proportion to be an 0{1) number but we are not aware 
of any reliable calculation of the ratio.) 


VI. GLUEBALL METASTABILITY AND RELIC 
DENSITY 

If the twin QCD sector had no couplings to other 
states, then after the phase transition, the glueball bath 
would behave roughly as non-relativistic strongly-self- 
interacting DM (since all of the glueballs have mass 
> 5.6Tc). That is, while number-changing (e.g. 3 — >■ 2) 
interactions were still fast enough to maintain number 
equilibrium, its temperature would decrease only log¬ 
arithmically with the scale factor |33) . with its energy 
density decreasing as p oc l/(a^ log(a/ao)), where oq is a 
constant. 

However, even in the absence of light twin sector states, 
the Higgs mixing portal with the SM necessarily provides 
a coupling to lighter states, and in particular, means 
that most of the glueballs decay rapidly to the SM. The 
question then becomes whether the relic population of 
(meta)stable glueballs is small enough not to be cosmo- 
logically dangerous. 

To calculate the relic population of (meta)stable glue¬ 
balls, we need to know the point at which their number¬ 
changing interactions freeze out. The last such interac¬ 
tions to freeze out will be processes with two particles 
in the initial state. Inspecting the glueball spectrum, if 
states have their equilibrium abundances then we expect 
the fastest such processes to be annihilation to lighter 
glueballs, e.g. 0 ^0 ^ —>■ We can therefore 

perform the usual calculation for the relic density of an¬ 
nihilating DM. If we parameterise the annihilation cross 
section as {av) ~ (^/(Aqqq)^, then the present-day relic 
density is (assuming no significant entropy injection) 

~ 3.8 X 10-i°(A(5CD/GeV)2C'-i , (6) 

in comparison to the present-day DM relic density 
— 0.12 [33]. This assumed that the SM and 
glueball temperatures were equal at the freeze-out time. 
This will be the case when the decay rates of some of the 
glueballs are sufficiently fast—roughly, faster than the 
Hubble rate at freeze-out. If it is not the case, then as 
described above, the glueballs will be at a higher temper¬ 
ature than the SM, as their temperature will have fallen 
only logarithmically with the scale factor. Comparing 
the —>■ SM decay rate (from the mixing of EqQ) 

to the Hubble rate at freeze-out, the former is larger for 


Aqcd ^ 0.6 GeV, so we are in the fast-decay regime for 
most of the Aq^^^ range of interest. 

As H/Hdm — few x I0 “®(Aqqj 3/GeV)^^”^, the relic 
density of stable glueballs will have no sigificant grav¬ 
itational effects, and if the metastable glueballs decay 
well before recombination time (~ 10^^ s), they will not 
inject enough energy to observably disrupt BBN or the 
GMB spectrum [35lI37j . However, an energy injection 
of > 10“^° of the DM energy density can, depending 
on the injection time and channels, have observational 
consequences if it occurs around recombination time or 
later, either through GMB effects, or via cosmic ray 
observations [371 |3H]. Thus, it may be a requirement 
that the meta-stable glueballs have lifetimes shorter than 
~ 1013 s. 

As an aside, note that the situation is different if 
mb' is light enough that the lightest twin QCD states 
are mesons rather than glueballs. The lightest meson 
state is a 0 so in the absence of lighter twin sector 
states, will decay through higher dimensional operators, 
as per the 0 glueball discussed below. The differ¬ 
ence is that its annihilations must now produce SM fi¬ 
nal states, rather than purely twin-QCD states, so will 
be much too suppressed to give it a sub-DM abundance. 
Number-changing interactions (e.g. 3 —?> 2 processes) will 
reduce its density, since other twin QCD states can de¬ 
cay to the SM, but since we are working in the regime 
of meson masses significantly larger than 1 GeV, the 
number-changing interactions generally freeze out before 
the abundance can be reduced to sub-DM levels [55] . 
Thus, it appears that, in the absence of twin sector CP 
violation and lighter twin sector states, the lightest b' 
meson must decay before BBN time in order to avoid 
dangerous energy injection. Since there are dimension-6 
operators that could lead to this decay, this condition is 
in principle easily satisfied. 

We now turn to the question of the lifetimes of the 
metastable 0 and glueballs mentioned in Sec¬ 
tion H. Let us first consider the case where there exist 
new SM-twin-sector interactions in the UV completion 
of the TH theory that allow the 0 ^ and 1^ glueballs to 
decay. For the 0 ^ glueball the lowest dimension effec¬ 
tive operators conserving total CP that allow the glue¬ 
ball to directly decay to the SM sector are of dimension 7, 
e.g. operators of the schematic form gqsg x iv{G'G')/M^ 
(here G’ are the SU{2>)' field strengths, while q stands 
in for SM fermions). These lead to a lifetime, Tq-+ ^ 
10“^^ s (M/5TeV)®(3GeV/AQQp)'^, irrelevantly short 
for astrophysical purposes (though interesting for dis¬ 
placed vertices at the LHG) unless the scale suppressing 
the operator is raised to M P 500 TeV. On the other 
hand, the leading effective operators for the I"' glue¬ 
ball are of dimension 10, e.g. operators of the schematic 
form x d^G'G'G'/M^. These lead to a lifetime 

Ti+- ~ 10 s (M/5TeV)^^(3GeV/AQ(2Q)^3^ which is of 
interest for indirect detection signals as discussed in Sec¬ 
tion VH. 

On the other hand, if only those interactions present 
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in the IR effective theory are considered then decays of 
the 1'' glueball compatible with conservation of both 
angular momentum and CP involve combinations of on- 
and off-shell glueballs, which then decay to the SM via 
the Higgs portal. For example, 1"' can decay to 0++ 
and an off-shell O'*"'*' / h, m a. I = 1 state. Since both C 
and P are violated in this process, and the only interac¬ 
tions in the twin sector that violate both C and P are 
SU{2y, the decay of the 1"' glueball needs to proceed 
through twin weak interactions involving both axial and 
vector currents. We can estimate its decay rate taking 
into account the leading heavy quark corrections to the 
vector and axial currents, and assuming that the decay 
to SM final states is via off-shell O'*"'*' / h mixing. Using 
the (dominant) one-loop b' quark correction to the vector 
and axial 311(2)' currents [toiHT] : 

(7) 


SJ ^ = 43/21^,2 

( 8 ) 

the decay rate of the 1+ glueball is estimated to be 
(ignoring dimensionless numbers and assuming ^ 
mo), 


ri+- oc 


'r;,(r 


19 4 4 ^4 


(9) 


and since r/i(mg) ^ TOq ~ mo we see that ri+- oc TOq^. 
Thus the decay rate of the 1"' glueball depends very 
strongly on the mass of the glueballs, or equivalently on 
the twin confinement scale. For example, for //?; = 3, 


/ mf,/ \12 y 10® s for Aqqq « 1.5 GeV 
^ \ 10® sfor A(5CD«3GeV 

( 10 ) 

i.e. changing Aq(^p by a factor of 2 leads to a 7-order of 
magnitude difference in the lifetime of the 1“' glueball! 
All that can be said is that, for reasonable values of the 
parameters of the model, it appears that the 1^ glueball 
may be a long lived state. Since it ultimately decays to 
SM states, even a tiny freeze-out density of 1"' glueballs 
may be dangerous, specifically if the 1"' lifetime is close 
to or after recombination time. We therefore from now 
on assume that, in the case where effects from higher 
dimensional operators can be neglected, Aq^q and mf,' 
are such that ri+~ <C tcMB — 10^® s, preferring larger 
values of Aq^^ or smaller values of . 

We finish by emphasising that the stability or duration 
of metastability of the 0 ^ and 1'' glueballs depends 
sensitively on whether the twin neutrinos are heavy or 
not, the mass scales mh> and Aq^q, whether twin CP is 
conserved, and finally the nature of the interactions in 


the UV completion between the twin and SM sectors in 
addition to the Higgs portal that might allow the glue¬ 
balls to more quickly decay. 


VII. INDIRECT DETECTION 

As well as leading to DM-SM scatterings in direct de¬ 
tection experiments, the Higgs mixing portal between the 
twin and SM sectors may result in SM energy injections 
from astrophysical DM-DM annihilations, leading to po¬ 
tentially observable indirect detection signals. For the 
twin DM candidates considered in this paper, annihila¬ 
tions proceed to lighter twin sector states. For t' and C 
DM, the dominant annihilation channel is to a b'b' pair 
via s-channel Z' exchange, due to the colour factor en¬ 
hancement of the final states available. Since 311(3)' is 
confining, the b'b' state will fragment into some number 
of twin glueballs and mesons — in the heavy-quark limit, 
we would expect dominantly glueballs. 

As previously discussed, if the C are light, then most of 
the glueballs will decay down to the lightest 0++ state by 
Cv' emission, with the 0'*"^ state then decaying via mix¬ 
ing with the SM Higgs. If there are no light hidden sec¬ 
tor states, then some of the glueballs may be metastable, 
while others will decay via the SM Higgs. In either case, 
we expect some combination of invisible decay products 
— C or stable glueballs — and of off-shell Higgs parti¬ 
cles, h*. The latter will have some distribution of mass 
squared values determined by the glueball masses and 
mass splittings. 

For Aqqq greater than a few GeV, most of these h* 
should be above the bb threshold. In that case, the SM in¬ 
jection products will, to a good approximation, be a spec¬ 
trum of bh with energies determined by the fragmenta¬ 
tion process. For DM of mass todm ^ 63 GeV as consid¬ 
ered here, the most sensitive probes of such astrophysical 
energy injection are cosmic ray antiprotons and gamma 
rays. The most up-to-date constraints on antiproton in¬ 
jection from dark matter annihilations come from the 
AMS-02 experiment [HI |35]. This sees an antiproton 
spectrum which is higher, at energies > 20 GeV, than 
expected in some models of astrophysical secondary pro¬ 
duction and propagation. If one takes these observations 
as disfavouring those models, and places constraints on 
the DM antiproton contribution assuming the secondary 
production models that better fit the AMS-02 data, one 
obtains constraints such as those of [43] , which finds that 
DM annihilating to bb with a thermal freeze-out cross sec¬ 
tion is disfavoured below todm — 100 GeV. If instead one 
includes the lower-secondary-production models in one’s 
uncertainty, the upper bound on the DM cross-section 
increases by almost an order of magnitude, and DM an¬ 
nihilating to bb is still viable in the entire mass range 
(though at the lower end, necessarily contributing signif¬ 
icantly to the observed AMS-02 spectrum). 

The best gamma-ray constraints on DM annihilation 
to hadrons, in this mass range, come from the FERMI 






LAT instrument. In particular, a recent FERMI analysis 
of dwarf spheroidal galaxies [44] places strong constraints 
on DM annihilating to bb with a thermal freeze-out cross 
section, disfavouring NFW-proflle DM in these galaxies 
with moM < 100 GeV. Taking into account uncertain¬ 
ties as to the DM profile in these galaxies, and/or us¬ 
ing the limits from the better-constrained diffuse Milky 
Way halo |3S], relaxes these constraints, disfavouring 
todm < 20 GeV. 

In our case, the uncertainty in the spectrum of h* en¬ 
ergies, and so bb energies produced, as well as the un¬ 
certainty regarding the invisible decay fraction, mean 
that we cannot make precise predictions for comparison 
with published limits. However, both the softening of 
the bb injection spectrum, and the invisible decay frac¬ 
tion, will act to decrease the strength of these constraints 
relative to straightforward bb injection, and in combina¬ 
tion with the astrophysical uncertainties, this means that 
these models are not ruled out by current data. However, 
especially at low tuning, and so low f /v, they predict in¬ 
direct detection signals potentially within the range of 
future observations. 

Particularly interesting from the point of view of fu¬ 
ture measurements are the tentative hints of a GeV-scale 
gamma-ray excess from the Galactic Gentre (GC) |^IT7] 
in FERMI data, which have now been corroborated 
by the FERMI collaboration itself [48|. This excess 
has a spectral and spatial profile compatible with ~ 
40 GeV DM annihilating to bb, with approximately ther¬ 
mal freeze-out cross section (or with ~ 10 GeV DM anni¬ 
hilating to leptons). While our DM masses are somewhat 
higher, the softer spectrum expected may plausibly give 
a good spectral fit. The higher mass and invisible decay 
fraction of our model would have some difficulty replicat¬ 
ing the approximately-thermal annihilation cross section 
derived in current analyses, but systematic astrophysical 
uncertainties are such that this may decrease with fu¬ 
ture observations. Alternatively, if Aqqq is low enough 
that a substantial fraction of the h* produced have mass- 
squared below the bb threshold, then injection to t+t“ 
and to lighter quarks may be important, likely lowering 
the required annihilation cross section (e.g. |49|). The 
potential compatibility with this existing signal provides 
additional motivation for better understanding the twin 
QGD fragmentation process. 

Another interesting possibility is that the fragmen¬ 
tation process produces an appreciable number of 
metastable glueballs, which may have a long enough life¬ 
time to travel for astrophysically-relevant distances be¬ 
fore decaying. For this effect to be significant in galactic 
DM observations, the lifetime of the glueballs would need 
to be on kiloparsec (~ 10 ^^ s) scales (just shorter than the 
lifetimes constrained by the GMB power spectrum). 

A different and more striking signal may occur for 
smaller lifetimes — if the lifetime is > the radius of 
the Sun, Rq ~ 2 s, but not much greater than the 
Earth-Sun distance, then metastable glueballs produced 
by annihilations of captured DM in the Sun can es¬ 


cape before decaying, and their decay products could 
be detected in charged cosmic rays or gamma rays |50^ 
155] . The solar capture rate of DM scattering spin- 
independently (and momentum-independently) from SM 
nuclei is, taking a fiducial mass of muM = 100 GeV, 
Cq ~ I0^° s~^ 153] — once enough DM has ac¬ 

cumulated in the Sun, an equilibrium will be reached in 
which the annihilation rate is half of the capture rate, 
bann = ^Cq. Since we have a precise prediction for the 
scattering rate off SM nuclei, for given DM mass, this 
determines the annihilation rate in the Sun (for annihi¬ 
lation cross sections of the magnitude we are consider¬ 
ing, the equilibrium rate is reached in much shorter than 
the lifetime of the Sun). Ref. [5T| considers the decay 
of metastable DM annihilation products into leptonic fi¬ 
nal states, finding that FERMI gamma-ray observations 
place an upper bound of Fann/dec ^ 5 x 10 ^® s“^, where 
= g-'Ro/T'-r _ g-AU/ 7 T jg fraction of events, for 
a mediator of proper lifetime r with boost factor 7 , in 
which the decay occurs between the Sun and the Earth. 
This is in the scenario where DM of mass 100 GeV an¬ 
nihilates inside the Sun to two mediator particles, each 
of which then decay to two leptons. In our case, only 
some fraction of the DM annihilation energy will go to 
metastable glueballs, and the gamma-ray spectrum pro¬ 
duced by each of these will generally be softer, so we 
would expect the constraints to be less severe. How¬ 
ever, for lifetimes ~ I s, but less than a few times 10^ s, 
there is some possibility that future observations could 
place limits on this scenario, or see a signal. Addition¬ 
ally, in the same way that constraints on an excess of 
high-energy electron/positrons from the Sun have been 
used to set limits on scenarios with leptonic decays of the 
mediators |S3J [SB] , a similar analysis applied to an¬ 
tiprotons may provide stronger constraints on the Higgs- 
portal models considered here. 


VIII. EQUILIBRATION OF SECTORS 

If the I/' are light, then we do not expect any of the 
glueballs to have cosmologically relevant lifetimes. Also, 
for most glueballs the decays involving ly' will occur at a 
faster rate than decays involving the SM. Since inter¬ 
actions with SM states are very suppressed, this raises 
the possibility that most of the entropy in the twin QGD 
sector may be transferred to the v', and remain there de¬ 
coupled from the SM. As discussed below, if the iz' are 
light, then they must be effectively massless as regards 
early-universe cosmology. Still, they could potentially 
give a large DR contribution, which would be in conflict 
with cosmological constraints on the extra number ATVeff 
of effective neutrino species. 

The minimum iz' energy density we might get would 
arise from the zz' being in equilibrium with the SM un¬ 
til after the twin QGD phase transition. In this case, 
since there were 5 * ~ 75 effective relativistic degrees 
of freedom in the SM bath, and SM neutrinos decou- 
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pie from the SM bath when 5 *,sm = 10.75, we have 

Ti, ~ = 1-971/' at late times, assuming 

that the ly' are light enough to behave as DR. Writ¬ 
ing the energy density in DR at around CMB times as 

Pdr = (3 ANee)pt,^sM (where p^^sM = | P'r is 

the naive SM neutrino energy density), the v' give a con¬ 
tribution to the effective number of neutrino species of 
AiVeff ~ (1.9)“^ = 0.075. If the t' are also light, the con¬ 
tribution is doubled. Non-instantaneous decoupling of 
the SM neutrinos gives a contribution of ANes = 0.046 
from the SM alone. Planck has measured ATVeff = 
0.15 ± 0.2 [34] (the ±0.2 being a simplified summary of 
a rather complex set of constraints: see [31] for details), 
so such a twin sector contribution is entirely compatible 
with present-day constraints. However, future observa¬ 
tions, including improved astrophysical determinations 
of the Hubble constant |31|, large-scale structure sur¬ 
veys and a future cosmic-variance CMB polarisation 
mission [58] , may be able to measure Ngg to an accuracy 
of CT(Aeff) 0.05, potentially testing the presence of a 
twin v' bath. 

Lowering the temperature at which the phase transi¬ 
tion occurs will lower the number of effective degrees of 
freedom in the SM bath, so will raise the contribution to 
ATVeff, up to ~ 0.1 for temperatures just above the SM 
QCD phase transition. Below the SM QCD phase tran¬ 
sition, the minimal contribution becomes ~ 0.5, which 
is observationally marginal. However, we will see below 
that we need to take ^2.5 GeV for sufficiently fast 

equilibration in any case, which is well above Aqcd- 

If the v' are too heavy to behave as DR until CMB 
times, they will constitute an extra (possibly hot) DM 
abundance. For greater than a few keV, and less 
than weak scale values, the v' will constitute a warm/cold 
DM abundance with far too high a relic density. If mu' is 
smaller than that, but still greater than 0(10eV), then 
they will constitute too large a warm/hot DM abundance 
(e.g. (SD). So, if rriyi < mhl‘2,, it is required to be less 
than a few eV, the same applying to t'. 

Thus, if the v' are light, then they are certainly much 
lighter than the glueball mass splittings. In this case, as 
mentioned above, all of the glueballs are unstable, either 
against decays involving v'v' or h*. The question of in¬ 
terest then becomes how much energy density from the 
twin QCD bath ends up in the eventually decoupled z/' 
and SM baths, and whether the v' represent an observa¬ 
tionally significant DR abundance. 

In the event that all of the entropy density of the twin 
QCD bath ends up in the v', we can estimate the result¬ 
ing contribution to AN^s as above. If mb' > 5 GeV, then 
during the first stages of its Boltzmann suppression, the 
b' coupling to the SM is large enough to equilibrate the 
SM and twin gluon sectors. Thus, the scenario resulting 
in maximum v' energy density is that in which the SM 
and twin QCD sectors decouple shortly after the b' go 
non-relativistic. In this case, the entropy density of the 
gluon bath, corresponding to 16 effective relativistic 
dof, stays within the twin sector. Since all of this gets 


transferred eventually to the zz', then if (as expected) they 
have a quasi-thermal spectrum, they will have a temper¬ 
ature ((16 ± 2 X |)/(2 X |))^^^ — 2.3 times the a~^ 
scaled SM temperature before the phase transition. So, 
they contribute AN^s — (2.3/1.9)^ ~ 2.0. 

A AiVeff contribution this large is incompatible with 
the existing cosmological observations. In fact, even if 
only the entropy of the twin QCD bath just about the 
critical temperature is transferred to the z^’, then the re¬ 
sulting AiVeff ~ 0.4 is in tension with observations. It is 
therefore a requirement that most of the entropy in the 
unconfined phase of the QCD bath is transferred to the 
SM, rather than to the v'. Logically, this could occur in 
a number of ways: 

• If twin QCD interactions with the v' and SM are 
too slow to transfer significant entropy to either in 
the first Hubble time, then the only process that 
can become important is glueball/meson decays. 
Then, the energy transferred to the v' and SM is 
set by the overall branching ratio of these decays. 

However, in our scenarios, SU{2y mediated pro¬ 
cesses are generally faster than Higgs portal ones, 
and we expect the branching ratio to favour v' over 
SM, so this does not solve the problem. 

• If twin QCD interactions with both the z/' and the 
SM are fast enough, then v' o QCD' -o- SM inter¬ 
actions may be fast enough to bring the z/' and SM 
sectors (almost) to equilibrium, and so to (almost) 
restore the AN^g = 0.075 situation derived at the 
beginning of this section. (Note that direct interac¬ 
tions between z^' and SM states are extremely sup¬ 
pressed, since diagrams must involve S'[/(2)' gauge 
bosons and multiple Higgs’s. The zz' may have a 
small interaction with the twin Higgs, giving them 
a < eV mass, but this coupling is again far too 
small to bring about equilibrium.) 

It is not trivial for interactions during the confined 
phase to be fast enough, since the glueball bath is non- 
relativistic, and its equilibrium number density falls very 
fast as the temperature drops. Energy transfer from the 
glueball bath to the SM needs to be fast enough that it is 
still operating until the SM and n' have reached almost 
the same temperature. 

In the confined phase, the fastest processes transfer¬ 
ring energy QCD' —>• SM are glueball and meson decays, 
for example, the 0++ —)■ /i* —)■ SM decay discussed pre¬ 
viously, and possibly SM-QCD' scatterings. Parameter- 
ising the effective SM decay rate of the glueball bath as 
PG = —rG->.SMPG + - ■ •; and solving the Boltzmann equa¬ 
tions numerically shows that Tg-j-sm ^ 10i7c is required 
to significantly reduce the v' abundance, where He is the 
Hubble rate at the start of the phase transition - see Fig¬ 
ure Since i7 oc T^, and the decay rates go as higher 
powers of A'q^^ (e.g. To-r-f oc A^^^q), increasing A'q^^, 
and so Tc, makes it easier to fulfil this condition, which 
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is satisfied for Aqqq > 2.5 GeV for our fiducial parame¬ 
ters. Since the decay and scattering rates of most of the 
glueball and meson states have not been calculated, the 
precise bound is uncertain. However, because these rates 
increase as much higher powers of the mass scale than 
the Hubble rate does, such uncertainty in the rates at a 
given energy has a small effect on the bound. 



FIG. 4. Relative energy density contributions, after 311(3)' 
phase transition, of SM (blue) and i/' (red) baths, for given 
glueball-SM energy transfer rate rG-»SM (compared to Hubble 
rate at start of phase transition He). This calculation assumes 
that the SM, QCD' and i/' baths are in thermal (chemical) 
equilibrium at T = Tc, which will not be the case if the cou¬ 
pling is weak enough (that is, the small-F part of this graph 
would have p^i even larger in a realistic scenario). The dashed 
green line shows the number equilibrium density of u'. 

The QCD' SM energy transfer from a mixture of 
confined and unconfined phases, at the critical tempera¬ 
ture, should (by detailed balance) be at most that from 
a wholly-confined-phase mixture (assuming that the con¬ 
fined phase loses most energy from processes with only 
twin sector particles in the initial state). Figure as¬ 
sumes that the rate is that of the confined phase de¬ 
cay approximation, from the critical temperature down¬ 
wards. So, either Aq^^ > 2.5 GeV, in which case we are 
driven almost to equilibrium, giving a AA^eff contribu¬ 
tion of ~ 0.075, or else a large fraction of the entropy 
of the twin QCD bath goes into the v', giving a AiVeff 
contribution of ~ 0.4 — 2.0 in conflict with observations. 


IX. SUMMARY AND CONCLUSIONS 

We have considered the status of DM (without a 
matter-antimatter asymmetry) in models based on the 
TH mechanism mm. focusing on its simplest imple¬ 
mentation, the Fraternal Twin Higgs [5], that contains 
just SU(3y and 311(2)' gauge interactions and only a 
third generation of fermions. While such models may 
feature a first-order 311(3)' phase transition at tempera¬ 
tures around or below DM freeze-out, this transition does 
not result in significant entropy production. 


The most natural DM candidate is the twin tau lepton 
t' , which requires masses m-r' > mhl2 in order to pro¬ 
vide the observed DM abundance. In particular, we find 
rui-' « 63 — 130 GeV for ratios f /v = 3 — 5, which im¬ 
plies a very mild 20 — 8% tuning. DM scattering off SM 
nuclei happens via Higgs exchange, and direct detection 
signatures lie below LUX current bounds [SD] but within 
LUX future sensitivity m- A very natural possibility 
within the Fraternal Twin Higgs scenario is that of mul¬ 
ticomponent DM. In particular, DM made of both t' and 
zz' species arises naturally as soon as the masses of both 
states are of similar size, though the phenomenology of 
this scenario is very similar to that of single-r' DM, mod¬ 
ulo the lack of any AfVeff contribution. In the case where 
rrir' + > mw' , gauge bosons become stable and 

may significantly contribute to the DM density, leading 
to a scenario where DM is made of three different species 
of twin particles (r', v' and W'^). We find, however, 
that this three-component scenario requires large values 
of //u, that result in fine-tuning between 5—1%, substan¬ 
tially worse than the single- or two-component cases. We 
briefly mention the possibility of twin baryon DM, where 
the DM particle would be a A' baryon made of three h' 
quarks. The efficient annihilation of h'b' pairs via twin 
strong interactions requires extremely large quark masses 
if the observed DM abundance is to be reproduced. In 
particular, we find that masses mb' > 1 TeV are required, 
which translates into f /v > 30 for yb' « 0.2, and there¬ 
fore a fine-tuning worse than 0.5%. 

Regarding possible indirect detection signatures, an¬ 
nihilation of DM particles happens mainly to b'b' pairs, 
which mostly results in glueball final states as a result of 
the fragmentation process in the twin QGD sector. Glue- 
balls will decay both into SM final states and into invis¬ 
ible twin sector states (either to in the case of light 
z/', or to (meta)stable twin glueballs). The former will 
consist mostly of bb pairs, giving broadly WIMP-like in¬ 
direct detection phenomenology, though with a presently 
non-calculable injection spectrum. The masses of our 
DM candidates are large enough that there is no incon¬ 
sistency with current data, though future observations 
should probe relevant regions of parameter space. De¬ 
cays of metastable glueballs may also have cosmological 
/ astrophysical consequences, though the strong depen¬ 
dence of these lifetimes on Aq^q, and on the UV physics, 
means that precise predictions are not possible. 

In the case where m^/ is small, a non-zero contribu¬ 
tion to the number of effective neutrinos, ANes, is a 
prediction of the theory. This contribution depends on 
the value of the twin confinement scale and we find that 
for Aqqq > 2.5 GeV, the i/' bath remains in equilib¬ 
rium with the SM bath after the twin QGD phase tran¬ 
sition and AVeff — AiVeff^M ~ 0.075, a value compatible 
with current bounds but within reach of future measure¬ 
ments. If the v' and SM baths fell out of equilibrium 
before the twin QCD transition, all the energy in twin 
glueballs is damped into the v' bath, in which case we 
would find ATVeff — ANes^sM ~ 2, strongly disfavoured 
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by current bounds. In the case of heavy twin neutrinos 
(toj,/ > mhl2), there are no light twin sector states, so 
no contribution to AA^eff- 

We remark that, independent of twin sector DM op¬ 
portunities, our analysis highlights cosmological and as- 
trophysical constraints on TH related models. The ANes 
bounds are independent of whether there is a twin sec¬ 
tor DM abundance, while the relic density calculations 
of Sections IV and |V| demonstrate that, for some mass 
parameters, the minimal TH model naively produces a 
super-DM abundance of stable states, which would need 
to be reduced by introducing new decay operators or hid¬ 
den sector states. 

Finally, we note that similar investigations of DM in 
Twin Higgs models have been carried out by other groups 

[SO] and ISIj. 
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